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ABSTRACT 


This  report  describes  exploratory  work  done  on  a  new  phase -shifting  technique. 
The  primary  objective  of  this  program  was  to  develop  and  demonstrate  a  laboratory- 
model  140-GHz  phase  shifter. 


The  desired  variation  in  phase  shift  is  produced  by  electronic  modulation  of  the 
width  of  a  rectangular  waveguide  which  results  in  a  change  in  the  phase  shift  per  unit 
length  along  the  waveguide.  The  change  in  effective  width  of  the  waveguide  is  accom¬ 
plished  by  means  of  a  PIN  diode  that  is  literally  distributed  along  a  sidewall  of  the 
waveguide. 


Among  the  tasks  requiring  study  were  1)  propagation  analysis  for  the  diode- 
loaded  waveguide  to  determine  the  relevant  parameters,  2)  fabricating  and  evaluating 
the  distributed  PIN  diodes,  and  3)  forming  a  phase  shifter  using  these  diodes.  Com¬ 
puted  and  measured  data  are  presented  for  waveguides  loaded  with  silicon  slabs  or 
various  conductivities  and  for  distributed  PIN  diodes. 

This  work  has  demonstrated  the  potential  feasibility  of  the  distributed  PIN  phase 
shifter. 
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Section  I 


INTRODUCTION 


The  primary  objective  of  this  program  was  to  develop  and  demonstrate  feasibil¬ 
ity  of  a  laboratory  model  140-GHz  phase  shifter.  The  significance  of  such  a  develop¬ 
ment  is  that  the  availability  in  some  practical  form  of  an  electronically  controllable 
phase  shifter  would  make  possible  the  construction  of  a  phased  array  type  of  antenna, 
thereby  providing  rapid  electronic  scanning  of  the  antenna  beam.  One  of  the  attrac¬ 
tive  features  of  millimeter  wavelength  radiation  is  that  it  can  be  focused  to  a  quite 
narrow  beamwidth  by  a  relatively  small  antenna.  The  high  resolution  (for  a  non- 
optieal  system)  and  small  size  of  such  an  antenna  point  to  applications  where  portable 
or  airborne  equipment  is  required  for  meteorology,  surveillance,  secure  communi¬ 
cations,  high-resolution  (by  microwave  standards)  imaging,  and  tracking  or  weapon 
delivery. 

The  approach  taken  here  to  millimeter  wave  phase  shifting  is  a  technique  in 
which  the  desired  v  riation  in  phase  shift  is  produced  by  electronic  modulation  of 
the  width  of  a  rectangular  waveguide,  which  results  in  a  change  in  the  phase  shift 
per  unit  length  along  the  waveguide.  The  change  in  effective  width  of  the  waveguide 
is  accomplished  by  means  of  a  PIN  diode  that  is  literally  distributed  along  a  sidewall 
of  the  waveguide  as  shown  in  Fig.  1.  With  no  bias  applied  to  the  diode,  the  transmis¬ 
sion  line  is  a  rectangular  waveguide  of  width  a,  loaded  with  a  dielectric  slab  of  width 
t.  Application  of  a  suitable  bias  signal  causes  injection  of  free  charge  carriers  into 
the  I  region  of  the  diode  in  sufficient  density  that  its  millimeter  wave  conductivity  is 
raised  to  a  high  value,  thereby  moving  the  effective  waveguide  wall  to  the  free  sur¬ 
face  of  the  diode  and  changing  the  width  of  the  waveguide  to  d  in  Fig.  1.  The  result 
is  a  change  in  the  waveguide  propagation  constant  and  therefore  a  change  in  phase 
shift  per  unit  length  relative  to  the  unbiased  condition. 

Among  the  problems  requiring  study  were  propagation  analysis  for  the  diode- 
loaded  waveguide  to  determine  the  rel<  vant  parameters,  fabrication  techniques  for 
use  and  evaluation  of  the  distributed  PIN  diodes,  and  forming  of  a  phase  shifter  using 
these  diodes. 

The  program  performance  goals  for  the  phase  shifter  were  as  follows: 

Frequency  140  GHz  Power  handling  capability  1  W,  CW 

Phase  shift  100° /cm  Drive  power  10  mW 

Insertion  loss  0.5  dB/cm  Modulation  rate  1  kHz 
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Although  development  did  not  progress  to  this  point,  no  fundamental  problems  were 
iound  in  the  proposed  phase  shifting  scheme.  Indeed,  the  technique  was  shown  to  be 
valid,  and  indications  arc  that  the  above  listed  numbers  are  realistie. 


Note:  The  diode  is  insulated  from  the  waveguide  walls  except  lor  its 
lower  terminal. 


Fig.  1.  Phase  shifter  geometry. 
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Section  n 


DEVICE  PROPAGATION  CHARACTERISTICS 

A.  LOADED  WAVEGUIDE  PROPAGATION  ANALYSIS 
1.  Introduction 


In  the  distributed  PIN  diode  phase  shifter,  the  diode  h:is  the  effect  of  a  semi- 
conductive  slab  loaded  in  a  section  of  rectangular  waveguide.  In  order  to  obtain  a 
good  understanding  of  the  device  operation,  a  theoretical  analysis  of  the  propagation 
in  loaded  waveguides  was  under uiKcn.  The  initial  results  of  this  analysis  have  been 
presented  in  an  RCA  internal  report.  In  this  section  we  present  an  extended  version 
of  that  analysis  and  remove  some  of  its  limitations. 


The  generalized  slab-loaded  waveguide  considered  here  is  shown  in  Fig.  2. 
Each  of  the  regions  is  assumed  to  be  uniform  and  homogeneous,  with  dielectric  con¬ 
stant  e j  and  conductivity  orj  .  The  previously  reported  analysis  is  a  special  case  of 
the  geometry  of  Fig.  2.  There,  region  1  was  the  semiconductor  slab  and  regions  2 
and  .‘5  were  free  space,  as  shown  in  Fig.  3.  In  that  analysis  the  phase  shift  of  the  de¬ 
vice  was  found  to  increase  considerably  with  increases  in  slab  thickness  up  to  a  cer¬ 
tain  critical  thickness.  Beyond  the  critical  thickness  the  analysis  essentially  broke 
down.  For  these  thiek  slabs  (thickness  greater  than  a  quarter  wavelength  in  tl  e  slab) 
it  appears  that  the  dominant  mode  electric  field  becomes  concentrated  in  the  semicon¬ 
ductor  poi'tion  of  the  waveguide  and  decays  exponentially  in  the  unloaded  section.  For 
slabs  under  a  quarter  wavelength  thick,  the  ''leetric  field  depend*  nee  is  sinusoidal 
both  in  the  loaded  and  unloaded  sections  of  the  waveguide.  The  present  analysis  yields 
the  propagation  constant  both  for  the  thinner  and  thicker  slabs.  However,  there  is  a 
definite  discontinuity  in  performance  at  a  thickness  greater  than  a  quarter  wavelength, 
with  the  thinner  slabs  yielding  improved  phase  shifi,/ insertion  los3  ratios. 

In  addition  to  treating  the  thicker  slab,  the  present  analysis  has  also  yielded 
results  for  other  slab  geometries.  This  includes  a  semiconductor  slab  mounted  off 
the  narrow  wall,  with  both  an  air  and  dielectric  backing.  Also  treated  is  the  ease  of 
two  semiconductor  slabs  symmetrically  mounted  on  the  two  narrow  walls  of  the 
waveguide. 


:s 


Fig.  3.  Single  slab  on  waveguide  sidewall. 


2.  Thcorctieal  Analysis 


In  studying  wave  propagation  in  hullow  conducting  waveguides,  it  is  common 
to  classify  the  solutions  into  two  types: 

1)  Waves  that  contain  eleetric  field  but  no  magnetic  field  in  the  direction  of 
propagation.  Since  the  magnetic  field  lies  in  the  transverse  plane,  the 
waves  are  referred  to  as  transverse  magnetic  (TM)  waves. 

2)  Waves  that  contain  magnetic  field  but  no  electrie  field  in  the  direction  of 
propagation.  Sinee  the  electrie  field  lies  in  the  transverse  plane,  these 
waves  are  referred  to  as  transverse  eleetric  (TE)  waves. 


The  above  is  not  the  only  way  in  whieh  the  possible  w.ive  solutions  may  be 
divided,  but  it  is  a  useful  way  in  that  any  general  field  distribution  excited  in  an  ideal 
waveguide  may  be  divided  into  a  number  (possibly  an  infinite  number)  of  the  above 
types  with  suitable  amplitudes  and  phases.  Of  eourse,  only  one  of  the  possibly  infi¬ 
nite  number  may  propagate  if  it  alone  is  excited  and  if  conditions  are  favorable  for 
its  propagation. 

The  normal  modes  of  propagation  for  a  waveguide  loaded  with  a  dielectric 
slab  are  not,  in  general,  TE  or  TM  modes,  but  a  combination  of  a  TE  and  a  TM  mode. 
However,  the  normal  modes  are  e  tner  transverse  electric  or  transverse  magnetic 
with  respect  to  the  slab  boundary  phones.  This  is  due  to  the  fact  that  normal  to  the 
slab  boundary  the  "guide"  appears  uniform,  with  the  narrow  walls  acting  as  a  short 
circuit  load.  The  field  components  are,  in  this  case,  derived  from  a  single  vector 
potential  which  is  directed  normal  to  the  slab  boundary  plane.  The  resultant  fields 
are  either  transverse  electric  or  transverse  magnetic  to  the  interface  normal.  For 
the  former  the  electric  field  lies  in  the  longitudinal  Interface  plane,  and  the  mode  is 
referred  to  as  a  longitudinal-seetion  electrie  (LSE)  mode.  For  the  latter  the  mag¬ 
netic  field  lies  in  the  longitudinal  interface  plane,  and  the  mode  is  referred  to  as  a 
longitudinal-seetion  magnetic  (LSM)  mode.  A  speeial  case  of  the  LSE  mode  is  the 
1  Eno  of  rectangular  waveguide.  These  modes  are  the  natural  modes  for  the  geom¬ 
etry  of  Figs.  2  and  3,  both  because  of  the  slab  orientation  and  also  because  the  fields 
will  be  launched  by  a  TE10  wave  in  empty  waveguide.  The  characteristic  equation 
for  the  propagation  constant  in  a  uniformly  slab-loaded  waveguide  is  transcendental, 
but  ean  be  solved  by  numerical  techniques. 

Propagation  equations  for  the  slab-loaded  waveguide  are  developed  by 
Collin  for  lossless  dieleetries.  These  expressions  can  be  modified  for  the  lossy 
case  by  replacing  the  (real)  dielectric  constant  by  a  complex  one  that  takes  the  con¬ 
ductivity  of  the  slab  into  aeeount.  For  the  LSE  mode  the  fields  ean  be  derived  from 
a  magnetic  Hertzian  potential  of  the  form 

n  =  (x,y)e_yz 
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fleKs  are  ‘°  ‘h°  S'ab-a,r  inlorf““-  The  and  magnetic 


E  =  jwaoVxn 


h  Vx  Vx  n 


Propagation  along  the  guide  is  given  above  by  the  exponential  factor  e~yZ  and  must 

mtions  ^h  m  * 1Cf1°aded  and  unloaded  portions  of  the  waveguide  if  the  boundary  con¬ 
ditions  at  the  interfaces  are  to  be  satisfied  for  all  arbitrary  values  of  z.  Here, 

Its  real  oart  ^  u  fPr°PaPtion  constant  composite  waveguide  structure, 

oh.  h  ’  h  attenuation  per  unit  length,  while  its  imaginary  part,  8,  is  the 

p  ase  shift  per  unit  length.  These  quantities,  actually  their  ratio  8  fa,  are  useful  in 
describing  the  phase  shifter  performance. 

wn),  .  5’onsider  a  sinK^  semiconductor  slab  mounted  on  the  narrow  waveguide 
a  as  s  own  in  Fig.  3.  Applying  the  boundary  conditions  that  the  fields  are  con- 

vanishes  o^th  1  °  lfnt°faoe  and  that  the  tan"ential  component  of  the  eleetric  field 
c  mshes  on  the  perfectly  conducting  guide  wa’ls  yields  the  following: 

h  tan  (it)  -  ~  l  tan(hd)  nv 


y2  =  2  -k  k  2  -  h2  -k  2 

O  o 


where 


2  9 

k  =  O)  u  € 
o  o  o 


k  ~  cr  “j  (cr/Wo) 


E5£:  3  and  4  h  and  £  are  the  transverse  (in  the  direction)  wave  numbers.  They 
are  the  two  unknowns  in  these  equations. 

,.Combining  Eqs*  3  and  4  yidds  the  following  transcendental  equation  in  the 
propagation  constant  y : 


0  =  f  (y)  =  ( y2  +  k  2 


i/2  [-/ 

tan  ( y2  +  k  k 


+  |y2  +  k  kQ2 


a  2 

tan  hy  +  k 


')  1 

•n 


(5) 


initially,  the  roots  of  Eq.  5  were  obtained  by  using  a  combined  Rayleigh-Ritz  Newton- 
Raphson  iteration  technique.  2  This  method  was  found  limited  in  its  application  and 
has  been  superseded  by  use  of  a  roct-seareh  algorithm.  The  computer  programs  for 
the  latter  approach  are  given  in  Appendix  A. 

In  the  root-search  method  of  solution  we  first  set  or  =  0  and  find  all  the  values 
of  and  0)  that  mal.a  f(y)  0.  Program  ROOT  of  the  Appendix  is  used  for  this  initial 
seareh.  Next  a  is  increased  from  zero  to  a  small  value,  0, 1  (ohm-m)"1,  and  the  ac- 
eurate  root  for  0  -  0  is  used  as  the  "first  guess"  for  the  root  of  Eq.  5  for  the  new  value 
°  '  The  new  root  is  determined  to  the  required  accuracy  by  incrementally  search¬ 

ing  the  a  -0  plane  in  the  neighborhood  of  the  "first  guess"  point.  This  process  is  eon- 
tinued  untR  a  reaches  a  high  value,  10 5  (ohm-m)-l,  and  then  repeated  for  various 
slab  thicknesses.  The  seareh  for  non-zero  values  of  a  is  accomplished  by  Program 
WAVE  in  the  Appendix. 

r  f u  The  root- seareh  method  used  in  the  computer  analysis  is  a  modified  version 

of  the  method  of  keenest  descent.  The  fact  that  y  is  a  complex  quantity  for  non-zero 
values  of  slab  conductivity  leads  to  a  two-dimensional  seareh.  The  technique  used 
was  to  step  both  m  and  0  in  both  the  positive  and  negative  directions  and  move  in  the 
direction  that  minimized  f(y)  in  Eq.  5.  When  the  current  step  size  did  not  lead  to  a 
smaller  value  of  the  function,  the  step  size  was  decreased  by  a  factor  of  ten.  This 
was  continued  for  five  iterations,  the  initial  step  size  in  a  and  0  being  1.0  and  the 
inal  step  size  being  0.0001.  In  the  computer  analysis  the  units  of  a  are  nepers/ 
meter  while  those  of  0  are  radians/meter.  Thus,  in  these  units,  the  values  of  ct  and 
p  are  accurate  to  the  fourth  decimal  place. 

For  a  =  0  the  system  is  lossless  and  the  propagation  eoi  stant  is  purely 
imaginary,  y  =  j0.  The  values  of  0  that  cause  f(j0)  of  Eq.  5  to  equal  zero  are  shown 
in  Fig.  4  as  a  function  of  relative  slab  thickness.  As  stated  previously  these  values 
ot  p  were  generated  by  Program  ROOT.  The  geometry  is  that  g*ven  in  Fig.  3,  at  a  fre¬ 
quency  of  140  GHz  in  standard  RG-138/U  waveguide.  The  branch  in  Fig.  4  marked 
dominant  mode  is  an  extension  of  the  TEjg  mode  in  unloaded  waveguide  (i.e,,  t/a=0). 

1  e  branch  marked  "0  =  k0"  is  a  trivial  root  of  f(j0).  However,  the  point  where 
tbe  dominant  mode  crosses  this  horizontal  line  corresponds  to  a  slab  thickness  of 
one-quarter  wavelength  in  the  dielectric.  It  will  be  shown  below  that  this  is  a  criti¬ 
cal  thickness  in  terms  of  phase  shifter  performance. 

.  F  ’ TS!n0,r  ,h0  lrlVia!  r00ts  kO’  F™m  Eq.  4  it  can  bo  soon  that  for  0<k„, 

Tho  »an  "'k  ,o  rca*; 'vhll°  for  ko  h  will  be  imaginary  anti  t  will  be  real 

Iho  wave  ratmbers  land  h  Eire  the'V  dependence  of  the  fields  in  the  semiconductor 
slab  and  air  portions  of  the  waveguide,  respectively.  A  real  wave  number  corresponds  to 
a  sinusoidal  dependence,  whereas  an  imaginary  wav0  number  corresponds  to  a  hyper- 
of  «Ponentlal  dependence.  Thus  below  the  "0=  Inline  the  dominant  mode  has 
sinusoidalfield  dependence  in  both  regions  of  the  guide;  while  above  this  line  the  sinusoidal 
ependonee  is  still  maintained  in  the  slab;  however,  the  fact  that  b  is  imaginary  ieads 
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Fig.  4.  Phase  constant  j3  forar  =  o  in  the  geometry  of  Fig.  3, 
as  a  function  of  slab  thickness. 


to  an  exponential  decay  i’i  the  air  region.  This  results  in  a  surface  type  wave  en¬ 
countered  in  microstrip  circuits  where  the  fields  are  concentrated  in  the  dielectric 
and  decay  exponentially  away  from  Its  surface.  3 


3.  Analytical  Results 


Using  the  various  roots  in  Fig.  4  as  a  starting  point,  Eq.  5  was  solved  to 
obtain  values  of  y  =  a  +  j0  for  the  case  of  non-zero  conductivity.  Initially  the  geom¬ 
etry  shown  in  Fig,  3  was  treated.  The  more  general  geometry  at  Fig.  2  will  be 
discussed  in  Sec.  II. A. 4. 

Typical  computer-generated  results  are  shown  in  Figs.  5  and  6  for  a  range 
of  slab  thicknesses.  Here  the  real  (attenuation)  and  imaginary  (phase  shift)  parts  of 
the  propagation  constant  are  plotted  as  a  function  of  the  slab  conductivity  with  slab 
thickness  as  a  parameter.  Initially  a  relative  dielectric  of  12  was  used;  later  results 
are  for  er  =  11.7,  which  is  more  ‘epresentative  of  silicon.  (This  change  does  not 
affect  the  results  significantly.)  In  these  curves  the  reference  for  phase  shift  is  taken 
to  be  the  zero  conductivity  slab.  The  insertion  phase  shift  for  a  section  of  waveguide 
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i  tjuito  largo  (.160  per  guide  wavelength,  about  0. 1  inch).  For  the  zero  eonductivity 
slab-loaded  waveguide  the  insertion  phase  shift  is  given  in  Fig.  4.* 

Examination  of  the  corresponding  curves  in  Figs.  5  and  6  leads  one  to  sus¬ 
pect  that  a  relationship  exists  between  a  and  0.  The  facts  that  the  real  part  peaks  at 
the  point  where  the  imagin  iry  part  goes  through  a  point  of  inflection,  and  that  the 
geneial  shapes  of  the  curves,  are  suggestive  of  the  transform  relationship  that  exists 
between  the  real  and  maginary  parts  of  various  physical  functions.4*5  However,  on 
closer  examination  of  the  curves,  it  appears  that  a  and  0  are  not  related  by  the  Hilbert 
transform  relationship  found  in  electrical  network  analysis.  This  appears  to  be  due 
to  the  fact  that  the  slab-loaded  waveguide  is  not  a  minimum  phase  shift  network. 

Further  effort  is  needed  to  determine  the  relationship,  if  any,  between  a  and  0. 

With  reference  to  i ig.  4,  the  results  in  both  Figs.  5  and  G  correspond  to 
initial  values  of  0  less  than  kQf  the  former  on  the  dominant  mode  braneh  and  the 
latter  on  the  lower  branch.  As  indieated  in  the  graphs,  the  eases  where  0  <  k0  yields 
useful  phase  shift/attenuation  characteristics,  while  for  0  >  kQ  the  attenuation  and  the 
phase  shift  increase  exponentially  with  slab  eonductivity.  These  results  agree  with 
the  above  postulation  that  the  millimeter  wave  fields  are  concentrated  mainly  in  the 
semiconductor  slab  when  0  >  kQ,  with  the  fields  being  highly  attenuated  as  the  slab 
becomes  more  lossy.  The  eleetrie  and  magnetic  field  distributions  in  a  dielectric 
(ferrite)  slab-loaded  waveguide  have  been  presented  by  Lax  and  Button.  ^  There  it  is 
also  shown  that  the  fields  become  concentrated  in  the  slab  as  its  thickness  is  increased. 

The  device  performance  was  evaluated  at  a  number  of  frequencies  in  the 
90-140  GHz  band.  Figure  7  shows  the  dominant  mode  phase  shift  and  attenuation 
for0<  kG  at  90  GHz,  whereas  Fig.  8  shows  the  maximum  (saturated)  value  of 
phase  shift  over  the  90  GHz  to  140  GHz  band  for  various  slab  thicknesses. 

All  of  the  analytical  resu  t.s  presented  above  have  been  for  standard  80-mil¬ 
wide  RG-138/U  waveguide.  In  Fig.  u  the  phase  shift  and  attenuation  are  plotted  for 
a  4.  8-mil-thiek  slab  for  two  different  waveguide  widths  (70  mils  and  90  mils).  These 
should  be  compared  to  the  t/a  =  G  percent  data  of  Fig.  5.  In  general  it  can  be  seen 
that  the  phase  shift  and  attenuation  increase  with  decrease  in  waveguide  width.  Thus, 
over  a  limited  range,  the  waveguide  width  provides  a  convenient  parameter  for 
adjusting  phase  shift  and  attenuation. 

1.  Analytical  Results  for  Additional  Slab  Geometries 

The  analytical  results  have  been  extended  to  a  number  of  interesting  slab 
geometries.  Among  them  are  a  semiconduetive  slab  mounted  off  the  narrow  wall 
of  the  waveguide  (with  and  without  a  dielectric  spacer  between  the  slab  and  the  wall) 


*The  conversion  from  radians/meter  to  degrees/cm  is  multiplication  by  a  factor  of 
0.573. 
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Fig.  9.  Phase  sliift  and  attenuation  for  single  slab  on  waveguide 
sidewall  at  140  GHz,  70  and  90  .nil  wide  guide. 

and  two  slabs  mounted  symmetrically  at  each  of  the  narrow  walls.  Both  of  these 
cases  can  be  treated  as  a  waveguide  containing  three  homogeneous  regions,  as 
shown  in  Fig.  2.  Application  of  the  boundary  conditions  yields  the  following 
transcendental  equation: 


0  =  f(y)  =  h^  tan  (h^)  t  tan  (h^)  +  h^  tan  (h^) 


-  h2  tan  (h1dj)  tan  (lyy  tan  (h2d2) 


(6) 


where 


and 


v2  . 2  2  .2  .  .  2  .2  .  ,2 
y  -h  -k  k  =  h  -  k  k  =  h  -k  k 

llo  2  2  o  23 


(7) 


\'Cr  t  =  1.2,3 

i 


(8) 
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The  method  of  solution  of  the  various  cases  is  similar  to  that  discusred 
above.  The  geometries  considered  here  were  not  too  different  han  that  shown  in 
Fig.  2.  This  similarity  allowed  the  use  of  the  roots  in  Fig.  4  as  a  "first  guess"  in 
the  iterative  solution.  The  search  technique  was  then  used  to  find  the  values  of  a  and 
0  that  eausc  the  right  side  of  Eq.  6  tu  go  to  zero.  For  this  case  a  variation  of  pro¬ 
gram  WAVE  is  used  in  which  Eq  7  is  substituted  for  Eq.  5. 

Typical  computer -gene rated  results  obtained  for  the  various  geometries 
are  shown  in  Figs.  10,  11  and  12.  For  the  eurves  in  Fig.  10  the  slab  width  is 
maintained  constant  and  the  slab  is  spaeed  progressively  off  the  narrow  wall.  These 
curves  arc  useful  for  tolerance  studies;  they  show  that  if  the  slab  is  located  slightly 
off  the  narrow  wall  the  deviee  performance  can  be  significantly  affected.  The  curves 
in  Fig.  11  correspond  to  a  semiconductor  mounted  on  a  dielectric  slab  (whose  dielec¬ 
tric  constant  is  the  same  as  that  of  the  semiconductor).  In  this  ease  the  sum  of  the 
two  slab  widths  has  been  maintained  constant.  These  curves  indicate  that  the  same 
amount  of  phase  shift  is  obtained  as  we  make  the  semiconductor  slab  thinner  but 
maintain  the  total  thiekness  cf  the  semiconductor  slab  and  dielectric  constant.  There 
is  a  limit  to  how  far  this  can  be  earried,  ir  that  the  semicond  letor  thickness  must  be 
at  least  one  skin  depth  in  its  high  conductivity  state.  Assuming  the  semiconductor  to 
be  quasi-metallie  in  the  high  conductivity  state,  the  skin  depth,  d,  ean  be  approximated 
by7 


_  _  (meters)  (9) 

"  f  1/2 

(n  f  ucr) 

where 

f  =  frequency 
H  =  slab  permeability 
0  =  slab  conductivity 

At  140  GHz  with  <j  =  104  (fl-m)“4  and  fi  =  u0,  the  corresponding  skin  depth  is  about 
0.053  mil,  or  about  1. 3  micrometers.  This  depth  corresponds  to  a  relative  slab 
thickness  of  less  than  0. 1  percent.  Another  requirement  on  the  composite  structure 
is  that  the  sum  of  the  semiconductor  and  dielectrie  slab  thicknesses  be  kept  less  than 
one-quarter  wavelength  in  order  to  avoid  the  surface-type  waves  discussed  above, 
where  the  field  becomes  concentrated  in  the  lossy  semiconductor.  One  possibility 
for  obtaining  the  two  layer  structure  is  to  epitaxially  grow  the  thin  semiconductor 
structure  on  a  high  resistivity  substrate.  The  epitaxial-silicon  technology  has  be¬ 
come  highly  advanced  over  the  past  decade,  and  appears  to  be  readily  applicable 
to  the  distributed  PIN  diode  phase  shifter  elements. 
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b.  Phase  shift. 

•  Phase  shift  and  attenuation  for  dieleetrie  spaeer  between  single 
slab  and  waveguide  sidewall  at  140  GHz. 


Tho  curves  in  Fig.  12  compare  the  symmetrically  and  unsymmetrieally  slab-loaded 
waveguides.  1  hese  curves  show  that  generally  splitting  the  slab  in  its  width  and 
mounting  the  two  slabs  symmetrically  on  opposite  narrow  walls  result  in  less  phase 
shift  as  compared  to  the  original  slab  mounted  unsymmetrieally.  The  lower  phase 
shift  results  from  the  fact  that  the  thinner  slabs  are  in  regions  of  weaker  eleetrie 
fierd  for  the  TE10  type  waves  that  arc  being  propagated.  The  symmetrical  arrange¬ 
ment  may  be  advantageous  when  using  two  semiconductor  slabs  whose  sum  is  greater 
than  a  quarter  wavelength  and  also  for  the  case  of  two  thinned  semiconductor  slabs 
mounted  on  lossless  dielectric  slabs.  For  this  ease  a  slab  whose  relative  thickness 
is  2  percent  of  that  of  the  'weguide  width,  mounted  on  14  percent  thick  dielectric, 
would  lead  to  surface  wave  i  in  the  semiconductor.  However,  making  the  arrangement 
symmetrical  (i.e.,two  1  percent  slab  on  7  percent  dielectrics)  might  lead  to  a  useful, 

efficient  device.  To  date,  these  structures  have  not  been  studied,  either  theoretically 
or  experimentally. 

Based  on  the  above  results,  the  optimum  geometry  (for  maximum  phase  shift 
and  minimum  drive  power)  at  140  GHz  in  RG-138/U  waveguide  is  a  semiconductor  slab 
of  1  percent  or  less  relative  thickness  mounted  or  grown  on  a  lossless  dielectric  whose 
relative  thickness  is  about  7  percent.  This  geometry  implies  a  semiconductor  slab 
thickness  of  less  than  0.8  mil  and  a  dielectric  thickness  of  5.6  mils.  This  structure 
would  yield  phase  shift  of  about  248° /cm  and  attenuation  due  to  the  semiconductor  of  a 
few  tenths  of  a  dB/em  in  both  its  low  and  high  conductivity  states.  Decreasing  the  wave¬ 
guide  width  would  lead  to  additional  phase  shift  with  increased  attenuation  as  shown  in 
fig.  9.  It  appears  th  t  two  of  these  composite  slabs  in  a  symmetrical  arrangement 
would  provide  about  500°  /cm  of  phase  shift  at  140  GHz  in  standard  RG-138/U  waveguide. 

5.  PIN  Diode  Geometry 


Two  sets  of  measurements  have  been  carried  out:  one  on  a  set  of  bulk  silicon 
slabs  of  various  conductivities,  the  other  on  a  set  of  distributed  PIN  diodes.  The  re¬ 
sults  of  the  bulk  slab  measurements  arc  reported  in  See.  II.  B,  while  the  PIN  diode 
results  are  given  in  See.  III.  For  the  bulk  slabs  the  geometry  approximated  by  the 
above  analysis  is  quite  good,  with  the  measured  results  for  both  thin  and  thick  slabs 
comparing  favorably  with  those  predicted  theoretically. 


However,  the  test  mount  used  to  evaluate  the  distributed  PIN  dio Jes  only 
approximates  the  structure  shown  in  Fig.  n-2,  in  that  the  diode  height  was  typically 
only  about  half  the  height  of  the  waveguide.  For  efficient  operation,  the  height  of 
the  diodes  is  limited  by  recombination  statistics.  The  excess  carrier  lifetimes  of 
the  distributed  PIN  diodes  has  been  running  about  3  us.  This  corresponds  to  an 
ambipolar  diffusion  length  of  3  mils.  (Diode  parameters  will  be  discussed  further 
in  See.  III.B.)  For  uniform  distribution  ot  the  injected  carriers  and  minimum  loss 
of  carriers  by  recombination,  the  width  of  the  diode  intrinsic  region  should  be,  at 
most,  a  few  diffusion  lengths.  This  limits  the  height  of  the  diodes,  for  the  geometry 
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of  Fig.  1,  to  about  10  mils.  To  partially  compensate  for  the  height  difference,  the 
waveguide  height  was  tapered  from  its  40-mil  standard  height  to  20  mils.  Any 
further  taper  was  found  to  lead  to  high  test  mount  insertion  loss  (>  10  dB)  due  to  the 
eoneentrated  nature  of  the  fields.  The  experimental  RF  result-  uve  discussed  in 
See.  III.C. 

The  ease  where  the  semiconductor  slab  only  partially  fills  the  waveguide 
height  is  difficult  to  analyze  in  that  it  leads  to  mixed  boundary  conditions.  However, 
based  on  the  experimental  PIN  diode  results,  it  appears  that  the  propagation  in  the 
half-height  loaded  guide  is  quite  different  from  that  in  full  height  loaded  waveguide. 
Thus,  the  latter  analysis  cannot  fully  explain  the  experimental  results.  A  possible 
analytical  approach  to  the  half-height  loaded  guide  is  to  break  up  the  waveguide  into 
a  number  of  uniformly  loaded  sections.8  Two  sueh  regions  are  needed  for  the  ease 
under  consideration,  as  indicated  in  Fig.  13.  The  solutions  for  the  fields  in  each 
region  are  obtained  (by  numerical  techniques,  if  neeessary)  and  the  continuity  of  the 
fields  across  the  boundary  of  each  region  yields  the  complete  solution.  While  this 
technique  is  rather  complex,  it  represents  a  systematic  approach  to  the  problem. 
Due  to  time  limitations,  this  analysis  has  not  been  carried  out. 

B.  BULK  SLAB  MEASUREMENTS 

Two  groups  of  bulk  silicon  slabs  were  prepared,  all  having  dimensions  40  mils 
high  and  1  em  long,  and  conductivities  covering  the  range  1  (ohm-m)-1  to  104 
(ohm-m)“4.  One  set  of  slabs  was  G.2  mils  wide  (t/a  =7.8  percent)  whereas  the 
other  set  was  13  mils  wide  (t/a  =  1G  pereent).  The  thicker  set  had  less  of  a  range 
in  conductivity,  thus  providing  limited  results.  The  slabs  were  placed,  one  at  a 
time,  along  the  sidewall  of  a  section  of  RG-138/U  waveguide  and  their  insertion 
loss  and  relative  phase  shift  were  measured. 

The  experimental  arrangements  used  to  make  these  measurements  arc  the 
same  as  those  shown  in  Fig.  24  for  PIN  diode  measurements,  except  that  the 
diode  bias  circuitry  is  absent;  attenuation  was  measured  by  the  substitution  tech¬ 
nique  while  phase  shift  measurements  were  obtained  by  means  of  a  standard  bridge 
technique.  Typical  results  obtained  for  the  two  sets  of  slabs  are  shown  in  Figs. 

14  and  15. 

The  insertion  loss  results  in  Fig.  14  are  in  good  agreement  with  the  similar 
data  of  Fig.  5  for  t/a  =  8  pereent,  whereas  the  measured  phase  shift  is  less  than 
predicted  by  the  analytical  study.  The  limited  experimental  data  for  the  1G  pereent 
slabs  in  Fig.  15  compare  qualitatively  with  the  "lower  branch"  analytical  results. 
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Fig.  13.  Geometry  for  regional  analysis  of  PIN  diode  loaded  waveguide 
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Fig.  14.  Measured  and  theoretical  phase  shift  and  insertion  loss  for  a 
7.  8 %  wide  silicon  slab  on  waveguide  sidewall  at  140  GHz. 
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15.  Measured  and  theoretical  phase  shift  and  insertion  loss  for  a 
16%  wide  silicon  slab  on  waveguide  sidewall  at  140  GHz. 
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Section  III 


THE  DISTRIBUTED  PIN  DIODE 


A.  FABRICATION 

I  he  distributed  PIN  diode  samples  used  in  this  investigation  were  formed  by 
means  of  conventional  diffusion  techniques.  The  major  aspects  of  the  fabrication 
process  are  outlined  in  the  following  paragraphs. 

1.  Intrinsic  Material 


The  "intrinsic"  material  from  which  sample  diodes  were  made  was  N-type 
silicon  having  various  resistivities  in  the  range  50  2000  ohm-em.  The  material 

was  obtained  from  commercial  sources  (Monsanto  and  Wacker)  in  the  form  of  wafers 
whose  thickness  would  constitute  the  height  of  the  diodes  from  contact  to  contact. 

The  range  of  thicknesses  was  G  to  13  mils.  The  silicon  was  of  good  quality, 
exhibiting  minority  carrier  lifetimes  in  the  range  100  to  1000ms  as  measured  on  the 
original  ingots  by  the  photoconduetive  decay  method.  (The  higher  lifetimes  correspond 
to  the  higher  resistivities. ) 

2.  Diffusion  of  P  and  N  Layers 


Iypieal  diffusion  cycles  were  begun  with  boron-doped  and  phosphorus-doped 
silicon  dioxide  deposited  with  silane  gas  at  350°C.  Capping  with  silicon  dioxide  was 
performed.  Diffusion  at  1200°C  for  3  to  0  hours  yielded  diffused  layers  10  to  15  (tm 
deep. 


An  alternate  diffusion  procedure  made  use  of  boron-doped  oxide  diffused 
at  lOSO^C  for  30  min.  This  boron-doped  Si02  was  left  on  to  act  as  a  cap  during  the 
phosphorus  diffusion.  The  diffusion  source  was  phosphorus  oxychloride,  at  a 
temperature  of  950°C  for  30  min.  This  procedure  yielded  P  and  N  layers  of  about 
2  to  3  Mm  thickness. 

Table  1  summarizes  the  details  of  the  diffusion  schedules  that  have  been 
used.  Additional  PIN  diffused  silicon  wafers  were  obtained  from  U.D.T.  Corp. , 
E.G.&  G.  Corp. ,  and  Unitrodc  Corp.  At  this  point  in  processing,  the  P,  I,  and  N 
layers  exist  over  the  whole  wafer,  but  individual  diodes  have  not  been  defined. 


TABLE  1.  DIFFUSION  SCHEDULES 


Schedule 

Diffusion  Conditions 

a 

1200°C  for  G  hours,  15-minute  quench  in  nitrogen. 

b 

1200°C  for  7  hours,  eool  to  800°C  in  4  hours,  reheat 
in  program -controlled  furnace  to  1200°  C  for  1  hour, 
eool  to  700°C  at  1  degree/minute,  quench. 

e 

1200°C  for  3  hours,  15-minute  quench  in  nitrogen; 
half  of  the  wafer  batch  slow-eooled  as  in  b.  above. 

d 

1050°C  for  1^2  hour  for  boron,  f>50°C  for  1/2  hour 
for  phosphorus. 

3.  Metallization 

Metal  electrodes  are  required  in  order  to  make  electrical  connection  to  the 
diodes.  These  electrodes  were  formed  by  sputtering  thin  (about  200°)  layers  of 
titanium  and  palladium  onto  both  sides  of  wafers,  followed  by  gold  plating  to  a  thick¬ 
ness  of  about  0. 1  mil.  The  gold  provides  the  bias  current -distributing  contact, 
while  the  titanium/palladium  layers  provide  adhesion  and  isolate  the  gold  from  the 
silicon  wafer.  This  isolation  is  necessary  because  gold  has  a  tendency  to  diffuse  into 
the  silicon  and  cause  a  substantial  reduction  in  carrier  lifetime.  For  the  same  reason, 
the  electrodes  must  be  removed  from  the  silicon  in  the  vicinity  of  the  paths  to  be 
followed  in  cutting  the  wafer  into  diodes,  lest  gold  be  driven  into  the  silicon  along  the 
edges  of  the  diodes.  This  last  step  is  accomplished  by  application  of  photoresist  to 
the  wafer  in  the  pattern  of  the  diodes  desired,  leaving  paths  several  mils  wide  between 
diodes.  The  electrodes  are  then  etched  away  from  these  paths,  leaving  the  wafer 
exposed  in  the  areas  where  cutting  or  dicing  is  to  be  done.  Figure  IGa  shows  a 
typical  wafer  having  patterns  for  two  sizes  of  diodes  on  its  face. 

4.  Dicing 

It  was  found  that  because  of  their  unusual  geometry,  typically  0.01  x  0.01  x 
0.4  in,  the  distributed  PIN  diodes  eould  not  be  separated  from  the  wafers  by  conventional 
scribing  techniques.  The  method  used  for  most  samples  was  to  cut  around  the  edges 
of  the  diodes,  to  a  depth  of  about  8  mils,  using  a  Quantronix  laser  seriber.  Figure  1Gb 
shows  a  view  of  the  diode  pattern  on  a  wafer  face  after  laser  scribing.  The  dark  lines 
are  the  laser-cut  channels,  and  the  wide  light-colored  strips  are  the  gold  contact 
strips.  Note  the  small  margin  of  silicon  around  the  contact  strips.  Since  the  laser 
seriber  was  limited  by  power  and  focusing  arrangements  to  a  cutting  depth  of  about 
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(a)  Wafer  with  patterns  applied 


(b)  Wafer  face  after  laser  scribing 


(c)  Finished  PIN  diodes  viewed  from  top  (contact) 
side;  three  widths  are  shown 


Fig.  10.  PIN  diodes  before  and  after  dicing. 


8  mils,  wafers  thicker  than  8  mils  required  that  the  iinal  separation  of  the  diodes  be 
accomplished  by  breaking  the  remaining  silicon  wf,b.  When  the  pattern  edges  and 
scribe  lines  followed  cleavage  planes,  separation  was  easy;  other  conditions  presented 
some  difficulty  leading  to  lower  yields.  The  final  step  in  processing  w'as  a  orief  acid 
etch  which  served  to  eliminate  undesired  surface  roughness  and  to  remove  surface 
imp  rities  picked  up  in  processing.  Figure  16c  shows  several  finished  diodes. 

An  alternate  method  of  dicing  that  has  been  tried  with  very  good  results  is 
based  on  the  use  of  a  multiple-wire  saw.  The  saw  that  was  tested  is  made  by 
Geoscience  Corp. ,  Mt.  Vernon,  N.  Y.  Using  a  fine  abrasive  slurry,  the  wire  saw 
gave  straight  smooth  cuts,  eliminating  the  problem  of  roughness  in  the  break  region  of 
laser  scribed  wafers. 

5.  Diode  Stacking 


Diode  stacking  refers  to  a  process  that  laminates  two  or  more  diode  wafers 
in  series,  thus  creating  a  monolithic  1  stack"  of  diodes.  This  technique  was  of 
interest  in  this  program  because  it  can  produce  an  overall  semiconductor  assembly 
whose  height  is  comparable  to  the  nominal  waveguide  height.  For  example,  a  stack 
of  5  diodes,  each  8  mils  thick,  will  produce  an  overall  diode  assembly  40  mils  high, 
thus  matching  the  nominal  waveguide  dimension.  Each  diode  would  operate  essentially 
as  described  for  the  individual  PIN  diodes.  Note  that  a  single  diode  with  a  very  large 
intrinsic  layer  will  not  perform  the  same  as  this  stacked  assembly.  Conductivity 
modulation  of  the  I  region  is  not  effective  much  beyond  the  diffusion  lengths,  or  5  to 
10  mils  in  practice. 

Diode  stacking  is  a  common  production  technique  for  high-voltage  rectifiers, 
and  the  type  of  diodes  being  investigated  for  the  phase  shifters  should  be  stackable. 

A  single  demonstration  effort  was  made  by  supplying  two  diffused  wafers  to  the 
appropriate  production  group.  About  500  A  of  platinum  were  sputtered  on  the  P  side 
of  one  wafer  and  the  N  side  of  the  other  to  be  joined.  Then  4000  4  of  titanium  were 
sputtered  over  the  platinum.  In  a  vacuum  fixture  the  plated  wafer  surfaces  were  then 
pressed  together  at  about  1000  psi  by  graphite  blocks  induction  heated  to  1200°C  for 
10  min.  Slow  cooling  at  l°/min  was  then  carried  out.  The  previously  described  wire 
sawing  technique  is  especially  applicable  for  dicing  these  laminated  wafers. 

B.  CHARACTERIZATION 

The  aim  of  this  section  is  to  characterize  the  PIN  diode  from  both  analytical  and 
experimental  viewpoints.  The  analytical  discussion  starts  with  the  continuity  equations 
for  electrons  and  holes  in  the  I  layer  of  the  diode  and  develops  the  current-voltage 
relationships.  The  experimental  discussion  reports  on  the  highlights  of  the  various 
measurements  performed  on  the  phase  shifter  PIN  diodes.  The  main  purpose  of  these 
measurements  is  to  characterize  diodes  that  will  provide  good  RF  performance. 
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Although  the  geometry  of  the  distributed  PIN  diode  is  somewhat  unique,  the 
conductivity  modulation  phenomenon  operates  in  the  same  manner  as  more  conventional 
PIN  diodes.  The  passage  of  moderate  amounts  of  forward  current  through  the  diode 
results  in  injection  if  holes  from  the  P-region  and  electrons  from  the  N-region  into 

e  central  intrinsic  region.  These  injected  carriers  significantly  increase  the 
conduct h  itv  of  the  I-region. 


'1  his  section  presents  those  aspeets  of  the  PIN  diode  theory  that  are  required  to 
unc  erst  and  the  operation  of  the  device  and  identify  the  significant  parameters  and  the 
ro  es  they  play.  The  purpose  of  the  following  paragraphs  is  to  analytically  determine 
he  cK odes  current-voltage  characteristic,  and  thereby  understand  the  details  of  the 
conductivity  modulation  phenomenon.  Essentially,  the  procedure  is  to  determine  the 
na  ure  of  the  charge  distribution  in  the  I  layer  at  the  various  injection- levels,  apply 
Poisson  s  equation  to  obtain  the  electrie  field,  and  then  integrate  to  obtain  the  voltage. 

Charge  Distribution  in  the  I  Layer 

i  -r  1  h?  StGady"State  ecIuations  relating  to  eurrent  flow  in  a  semiconductor  are 
the  drift  and  diffusion  current  equations  and  the  continuity  equations  for  electrons  and 

into'' the  following^6  diffUSi°n  dominated  Ration,  these  equations  can  be  manipulated 


Unr  dx2 


(9) 


J  -  e/j^  [(b+l)nE  +  5(b-l)dn/dx] 


(10) 


where  a  one-dimensional  model  has  been  assumed  without  saerifieing  the  generality 
of  the  results,  and  where 


i? 

k 

T 

e 

T 


=  the  equal  electron  and  hole  densities 

kT 

e 

-  Boltzmann's  constant 
=  temperature,  K 
=  electronic  charge 
=  high  level  lifetime 
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i  (oru  ) 
n  P 
b 

x 

J 

E 


electron  (or  hole)  mobility 

U  n4i 

n  p 

distance  into  I  layer;  x  =  o  at  P-I  junction  and  x  =  L  at 
N-I  junction 

current  density 

electric  field  intensity 


Equation  9  has  the  solution 

nQ  sinh  [(L-x)/L  ]  +  n  sinh  (x/L  ) 

n=  sinh  (L/L  ) - —  (11) 

a 


where 

1  /2 

La  =  t2^Mnr /(b+1)]  =  ambipolar  diffusion  length 

L  =  length  of  I  layer 

nQ  =  carrier  density  at  P-I  junction  (x=o) 
nL  =  carrier  density  at  N-I  junction  (x=L) 


at  thP  J  H  ™  6  ,he  constants  no  a"d  nL  we  apply  boundary  conditions 

two  junctions.  This  is  done  by  equating  currents  on  each  side  of  each  of  the 

junctions.  For  the  P-I  junction  it  is  convenient  to  consider  the  electron  current 
(  ensity)  Jn-  'The  low  resistivity  of  the  P  layer  results  in  only  a  small  electric  field 
which  can  be  neglected  leaving  Jn  as  a  purely  diffusive  current.  In  the  I  region,  J 

as  both  diffusion  and  drift  components.  Similar  statements  apply  for  the  hole  currents 
on  each  side  of  the  N-I  junction. 


For  the  P-I  junction  we  obtain 


0n 

o 

L  p 
np 


J  +  2  8  dn 

e/ip  (b+1)  (b+1)  dx 


x  =  o 


(12) 
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For  tbs  N-I  junction  we  obtain 


0  ni 


L  n 
P  n 


eu  (b+1) 
P 


2j9b  dn 
(b+1)  dx 


x  =  L 


(13) 


wh  re 

J  =  total  current  density 

L  (or  L  )  =  electron  (or  hole)  diffusion  length  in  the  P  (or  N)  layer* 
n  p 

p  (or  n  )  =  hole  (or  electron)  density  in  the  P  (or  N)  layer. 

P  n 

The  next  step  might  be  to  substitute  the  expression  for  n  of  Eq.  11  into 
Eqs.  12  and  13  and  solve  them  simultaneously  for  nQ  and  n^,  a  difficult  task.  How¬ 
ever,  two  approximations  are  applicable,  either  of  which  will  simplify  the  solution 
of  these  equations.  Solutions  for  nQ  and  nL  are  readily  obtained  for  L>>  La 
("long  diode")  or  L  <<  La  ("short  diode").  The  solutions  have  the  same  form  for 
either  case,  but  different  coefficients.  The  long  diode  approximation  allows  Eq.  11 
to  be  written  as 


n  =  n  exp  (-x/L  ) 
o  r  a 


(14) 


While  the  short  diode  approximation  gives 


n  =  n  "T(l“T)  +  n  — — 
o  2  '  W  L2L 


(15) 


Under  these  conditions  Eq.  12  and  13  become,  respectively,  for  the  long 


diode 


0n 


2$n 


L  p  eu  (b+1)  (b+1)  L 

n  p  P  a 


(16) 


*If  Ln  or  Lp  is  much  greater  than  the  length  of  the  relevant  layer,  the  length  of  the 
layer  replaces  the  diffusion  length  in  Eqs.  12  and  13. 
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J 


(17) 


2j3bn 


yv  eyb+i>  <b+i>  La 


and  for  the  short  diode 


3"°  J  8W 

LnPp  e4p(b+l)  (b+l)  L 


0", 


Lp"n  e«p(b  +  D 


B  b  <°r  -  V 

(b  +  1)  L 


(18) 


(19) 


If  we  now  consider  asymptotic  solutions  of  the  above  four  equations,  we  find  two 
solutions  for  each  equation,  corresponding  to  two  distinct  injection  levels.  For 
example,  in  Eq.  1G  we  can  consider 


or 


n  << 
o 


2  L  p 


(b+l)  L 


a 


("low  injection  level") 


n  >> 
o 


2  L  p 
n  JP_ 


(b+l)  L 


a 


("high  injection  level") 


Similar  relations  apply  for  Eqs.  17,  18  and  19.  Solving  the  four  equations  (16,  17, 
18,  19)  under  these  approximations,  we  find  for  both  the  long  and  short  diodes  that 
the  low  injection  levels  result  in  nQ  and  nL  proportional  to  J  while  the  high  injection 
levels  result  in  nQ  and  n£,  proportional  to  J^/2.  if  these  results  for  nQ  and  nj_,  are 
substituted  into  Eq.  11  we  obtain  the  I  layer  charge  distribution  for  low  and  high 
injection  levels  in  long  and  short  diodes. 

2.  I-V  Characteristic 


The  electric  field  in  the  I  layer  is  found  for  the  various  cases  by  substituting 
the  expressions  for  the  charge  distribution,  n,  into  Eq.  10.  For  long  and  short  diodes 
at  low  injection  levels  this  step  eliminates  J  from  Eq.  10  indicating  E  is  independent 
of  J .  The  voltage  drop  across  the  I  (bulk)  layer  is  found  by  integrating  E.  For  both 
diodes  the  voltage  is  constant,  given  by 
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(20) 


g  hl/2w 

V_  = — .  ■-  exp(L/2L  )  where  L>>  L 

B  (b+1)  a'  a 


V_  =  (L/L  )2  where  L<<  L 

B  (b+i,2 


(21) 


Since  the  injected  carrier  density  was  found  to  be  linearly  proportional  to  the  total 
current  during  low  injection,  the  condui  livity  must  be  inversely  proportional  to  the 
current  and  the  I  layer  voltage  drop  should,  indeed,  be  a  constant. 


For  high  injection  levels  the  injected  carrier  density  was  four.J  to  be  propor¬ 
tional  to  jl/2  Inserting  the  specific  relations  into  Eq.  10  we  find  that  E  is  no  longer 
independent  of  J.  Integrating  E  to  obtain  the  I  layer  voltage  we  find,  for  both  the 
long  and  short  diodes, 


J  -  A  VB2  (22) 

where  the  constant  A  has  appropriate  values  for  the  long  and  short  diodes. 

Figure  17  shows  the  two  injection  regimes  discussed  here.  In  addition, 
at  very  low  drive  levels  where  the  injected  carrier  density  is  less  than  the  intrinsic 
concentration,  the  I  layer  is  ohmic  in  behavior,  presenting  a  third  regime. 

This  analysis  deals  with  VB,  the  voltage  drop  across  the  I  layer.  To  accurately 
predict  the  total  diode  voltage  drop  we  must  add  to  VB  the  voltage  drops  across  the 
contacts,  the  doped  layers,  and  the  junctions.  While  the  analysis  given  here  is 
adequate  for  our  purposes,  its  approximate  nature  must  be  kept  in  mind. 

3.  Experimental  Measurements 


When  diodes  have  been  constructed  there  are  three  measurements  which  can 
be  made  with  reasonable  ease  to  help  determine  whether  the  diodes  will  function 
properly.  They  are  spreading  resistance,  I-V  characteristics,  and  carrier  lifetime 
measurements.  The  fabrication  details  referred  to  are  described  in  part  A  of  this 
section. 


a.  Spreading  Resistance 


Spreading  resistance  measurements  were  performed  on  the  various 
diffused  silicon  PIN  diode  wafers  by  using  a  four-point  probe  assembly.  After  the 
wafers  are  lapped  down  to  a  five-to-one  taper,  the  local  resistivity  is  measure:4  by 
using  the  probe  assembly.  These  results  are  summarized  in  Table  2.  Group  3  diodes 
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Fig.  17.  Theoretical  I-VB  characteristic  for  PIN  diode. 


TABLE  2.  SPREADING  RESISTANCE  MEASUREMENTS 


Group 

Diffusion 

Schedule 

L  ,  L 
n  P 

(U  m) 

P  *  n 

P  n 

(cm'3) 

L 

(Mm) 

P 

(  Cl -cm) 

3 

a 

13.5 

4  x  1019 

212 

50 

4 

b 

15.6 

5  x  1019 

208 

650 

9 

d 

2.5 

00 

rH 

O 

236 

1600 

NOTES 


L 

P 


are  the  lengths  of  the  doped  layers. 


p  ,  n  are  the  carrier  densities  of  these  layers, 
p  n  J 

L  is  the  length  of  the  intrinsic  layer. 

P  is  the  resistivity  of  the  intrinsic  layer. 
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were  processed  from  50  ohm-cm  silicon  wafers,  while  groups  4  and  9  wafers  had 
initial  resistivities  of  750  ohm-em  and  near  2000  ohm-em,  respectively.  Groups  3 
and  4  are  seen  to  have  fairly  similar  doped  layers,  while  group  9  diodes  have  a  more 
shallow  diffusion. 

The  spreading  resistance  data  has  been  useful  in  evaluating  the  post¬ 
diffusion  wafers,  by  providing  the  carrier  populations  and  lengths  of  the  various  regions 
of  the  PIN  diode  phase  shifting  elements.  It  has  confirmed  that  the  manufacturing 
processes  are  providing  diffusions  of  the  type  desired. 

b.  I-V  Characteristics 


The  different  diode  groups  exhibited  a  large  range  of  I-V  characteristics. 
Figure  18  is  a  multiple  exposure  photograph  of  the  I-V  characteristics  of  a  number  of 
diodes  from  various  groups.  Diodes  from  groups  3  and  9  are  seen  to  have  a  sharp 
break  in  their  characteristic  while  the  group  4  diode  is  seen  to  have  a  gradual  break 
and  a  corresponding  large  forward  voltage  drop.  Generally  diodes  with  a  sharp, 
nearly  vertical  break  at  about  0.6  V  gave  satisfactory  insertion  loss  and  phase  shift 
performance,  while  diodes  with  a  large  forward  voltage  drop  exhibited  high  insertion 
loss. 


G  90 2 

G  3D4— ^  J  G4DI 


V  E  R  T.  = 
50  mA/div 


HOR.  =  0.5  V/div 


Fig.  18.  Multiple  exposure  photograph  of  the  I-V 

characteristic  of  various  distributed  PIN  diodes. 


A  typical  plot  of  the  I-V  characteristic  of  a  group  9  diode  on  a  log  I-log 
V  format  is  shown  in  Fig.  19.  This  eurve  (which  plots  the  total  voltage  across  the 
diode)  is  to  be  compared  with  the  theoretical  eurve  shown  in  Fig.  17  (which  plots  only 
the  voltage  across  the  intrinsic  layer).  To  get  the  intrinsic  layer  voltage  from  the 
total  voltage,  we  would  have  to  subtract  the  junelion,  eontact  and  series  resistance 
voltage  drops.  However,  it  is  possible  to  evaluate  the  experimental  data  without 
actually  doing  this.  In  Figs.  17  and  19  there  are  a  number  of  essential  similarities. 
For  instance  at  low  values  of  current  and  voltage,  the  curve  in  Fig.  19  exhibits  a 
linear  region  (slope  <-a  1)  corresponding  to  the  ohmic  regime  of  Fig.  17.  This  region 
is  followed  by  a  steeply  rising  portion  corresponding  to  the  vertical  low  diffusion 
regime  of  the  theoretical  cui’ve.  The  falloff  of  the  steep  portion  in  Fig.  19  at  high 
values  of  diode  current  corresponds  to  the  high  diffusion  regime.  The  break  between 
these  two  regimes  is  seen  to  oeeur  at  a  few  hundred  milliamperes  of  forward  eurrent. 


Additional  insight  ean  be  gained  by  plotting  the  data  of  Fig.  19  on  a  log 
I-V  (i.e.,  semilog)  format,  as  shown  in  Fig.  29.  This  eurve  shows  a  straight  line 


eV 

exponential  regime  where  I  ~exp  (  jr^)»  where  n 


2.5  for  this  particular  diode. 


The  theoretical  value  of  n  for  a  PIN  diode  in  the  double-injection  regime  is  two. 

The  falloff  of  the  eurve  at  high  currents  characterizes  the  effeets  of  series  resistance. 
In  this  region,  the  eurrent  departs  significantly  from  exponential  behavior  and  becomes 
limited  by  the  series  and  eontaet  resistances  of  the  diode. 


II  is  possible  to  use  the  eurve  in  Fig.  20  to  get  a  measure  of  the  degree 
of  conductivity  modulation  under  bias  conditions.  The  diode  series  resistance  ean  be 
effectively  subtracted  by  extending  the  straight  line  section  of  the  eurve.  The  ae 
conductance  of  the  intrinsic  layer  is  then  given  by  the  slope  of  the  I-V  eurve  at  the  bias 
eurrent  of  interest.  Figure  21  shows  the  conductivity  modulation  as  a  function  of 
bias  current.  The  initial  conductivity  of  the  I  region  for  this  diode  was  about  0.  00 
(fl~m)”l.  Thus  we  have  effected  about  a  four-order-of-magnitude  modulation. 


e.  Exeess  Carrier  Lifetime 


An  important  parameter  in  characterizing  a  PIN  diode  is  the  lifetime  T 
or  the  ambipolar  diffusion  length  La  of  the  exeess  carriers  injeeted  into  the  I  region. 
These  two  parameters  are  related  by 

1  /? 

L  =  \2$n  T/(b+l)]  (23) 

a  n 

For  intrinsic  silicon  this  bteomes 

1/2 

L  (mils)  1.7  [T  (Ms)] 
n 


(24) 
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Fig.  21.  Conductivity  of  intrinsic  region  of  a  PIN  diode  as 
a  function  of  forward  eurrent. 

A  plot  of  Eq.  24  is  given  in  Fig.  22.  From  Eq.  11  it  ean  be  seen  that  La  is  the 
characteristic  length  for  PIN  diodes.  For  a  long  diode  (L»  La)  the  execss  carriers 
decay  as  exp  (-x/La)  from  the  injeeting  eontaet. 

The  cireuit  used  for  lifetime  measurements  of  the  PIN  diode  phase 
shifter  elements  is  shown  in  Fig.  22.  The  technique  used  is  to  fill  the  intrinsic  layci 
of  the  diode  with  cxeess  holes  and  electrons  by  pulsing  the  diode  under  test  (DU1) 
with  a  forward  voltage,  and  then  observe  the  open-circuit  voltage  transient  across  the 
diode.  Under  conditions  of  high-level  injection  (i.c.,  injected  carrier  density  greater 
than  the  doping  density  of  the  I  region),  the  lifetime  is  related  to  the  voltage  deeay  as 
follows 

-1 

2KT  dV 

T'—  (ir» 

Table  3  summarizes  the  lifetime  measurements  for  the  various  groups  of 
PIN  diodes.  Sinec  the  lifetime  is  a  function  of  the  injection  level,  it  is  important  that 
it  is  measured  at  the  level  of  interest.  The  values  reported  in  the  table  were  measured 
at  a  forward  current  of  300  mA,  and  range  from  about  2.5  Ms  to  4.  0  fJs.  From  Fig. 

22  these  values  correspond  to  a  range  of  ambipolar  diffusion  lengths  between  2.0  mils 
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Fig.  22.  Ambipolar  diffusion  length  versus  excess 
carrier  lifetime  for  intrinsic  silicon. 


and  3.4  mils.  Most  PIN  diodes  evaluated  have  approximately  a  10-mil  intrinsic  layer 
width,  vvhieh  corresponds  to  about  three  to  four  diffusion  lengths. 

An  inerease  in  T  by  about  one  order  of  magnitude  would  lead  to  a  9-mil 
diffusion  length,  which  would  increase  the  injection  efficiency  and  would  lead  to 
lower  drive  power  requirements. 

C.  II F  MEASUREMENTS 

1.  Test  Setup 

'1  he  experimental  arrangements  used  to  measure  attenuation  and  phase  shift 
are  shown  in  Fig.  24.  Insertion  loss  was  measured  by  the  substitution  method  by  using 
the  calibrated  variable  attenuator.  Plots  of  loss  versus  bias  current  were  obtained 
by  means  of  the  sweeping  current  supply  and  the  X-Y  recorder.  Use  of  the  log 
converter  allowed  calibration  of  the  recorder  charts  in  dB/in  on  the  Y  axis,  while  the 
X  axis  units  were  mA/in. 

Phase  shift  measurements  were  done  by  means  of  a  standard  bridge 
technique,  Fig.  24b.  Here  the  signal  from  the  klystron  source  is  split  between  two 
paths.  One  path  provides  the  phase  reference  signal,  while  the  other  includes  the 
phase  shifter  producing  unknown  phase  shift.  The  probe  of  the  slotted  waveguide 
section  is  moved  along  the  guide  until  a  null  is  encountered  in  the  standing  wave 
produced  by  the  signals  from  the  two  bridge  arms.  (The  variable  attenuators  may  be 
adjusted  for  a  deep  null.)  Next  the  phase  shifter  bias  current  is  increased  from  0  to 
500  mA  in  convenient  increments  while  the  location  of  the  null  is  tracked  with  the 
slotted  line.  Provided  the  reference  arm  has  not  been  disturbed,  the  difference  in 
phase  shift  between  any  two  bias  current  levels  can  be  obtained  in  degrees  by  dividing 
the  linear  null  shift  by  the  guide  wavelength  and  multiplying  by  3G0.  In  addition  to  the 
bridge  configuration  shown  in  Fig,  24b,  occasional  tests  were  done  by  using  an 
isolator  in  eaeh  arm  of  the  bridge  to  cheek  for  the  absence  of  reflection  effects. 

2.  Diode  Mount 


Because  of  the  distributed  parameter  nature  of  the  PIN  diodes  used  here,  the 
experimental  diode  mounts  were  essentially  nothing  more  than  appropriate  sections  of 
transmission  line.  Of  eourse,  the  practical  matters  of  matching  to  standard  waveguide, 
mounting  the  diode,  providing  a  bias  eurrent  path,  and  allowing  for  easy  installation 
and  removal  of  the  various  diodes  required  the  addition  of  appropriate  features. 

The  test  equipment  and  components  available  for  immediate  use  were  mostly 
in  RG138/U  waveguide,  which  has  an  interior  width  of  80  mils  and  a  height  of  40  mils. 
Because  of  diodes  were  limited  to  heights  in  the  5  to  15  mil  range  by  the  nature  of 
the  conductivity  modulation  mechanism,  the  height  of  the  waveguide  was  redueed  in 


the  diode  mount  to  20  mils  or  less  while  retaining  the  80  mil  width  for  convenience. 

The  diodes  were  held  in  place  in  the  mount  by  low  density  plyfoam  spacers  or  wedges 
of  Rexolitc  1422  inserted  in  the  waveguide.  The  dc  bias  circuit  was  established  by 
allowing  one  contact  strip  of  the  diode  to  rest  against  the  bottom  of  the  waveguide 
while  a  fine  insulated  wire  was  inserted  through  a  small  hole  in  the  top  of  the  waveguide 
to  make  contact  with  the  other  diode  contact  strip.  Figure  25  shows  a  cross-section 
view  of  the  diode  mount,  looking  through  the  narrow  wall  of  the  waveguide,  while 
Fig.  25a  is  a  photograph  of  the  same  mount. 

Several  mounts  were  m  .dc  and  tested  for  transmission  loss  without  diodes  at 
140  GHz.  The  2  to  3  dB  typical  measured  loss  is  not  considered  high  since  the  mount 
is  in  the  form  of  a  split  waveguide  embedded  in  a  brass  block,  and  has  an  overall 
length  of  4  in,  including  the  region  v/here  the  diode  is  located,  the  input  and  output 
tapers,  and  input  and  output  waveguide  extensions  for  working  convenience.  While 
the  2  to  3  dB  mount  loss  would  be  unacceptable  in  a  phase  shifter  for  actual  antenna 
use,  it  is  not  a  cause  for  concern  when  evaluating  the  PIN  diodes. 

3.  Test  Data 


Various  distributed  PIN  diode  samples  prepared  as  described  in  Sec.  III.  A 
were  tested  for  attenuation  and  phase  shift  at  140  GHz.  All  diodes  were  1  cm  long 
and  were  tested  for  forward  bias  currents  in  the  0  to  500  mA  range. 

It  is  important  to  note  that  almost  every  diode  tested  differed  in  some 
significant  way  from  the  model  used  for  theoretical  studies.  These  differences  were 
the  result  of  expediencies  of  the  moment;  a  practical  approximation  to  the  analytically 
studied  configuration  does  appear  to  be  attainable. 

The  test  results  wore  quite  sensitive  to  the  position  of  the  diode  within  the 
mount.  For  any  given  diode,  the  strongest  effects  were  due  to  the  proximity  of  the 
diode  to  the  waveguide  sidewall.  Figure  2G  shows  how  the  loss  characteristic  of 
diode  G19-D2  varied  for  1-mil  changes  in  the  diode  position  with  respect  to  the  side- 
wall.  The  indicated  spacings  were  obtained  by  means  of  IVTylar  shims  placed  between 
the  diode  and  the  waveguide  wall.  Figure  27  shows  similar  data  for  seven  different 
positions  of  diode  G19-1F. 

The  relative  heights  of  the  diode  and  the  waveguide  also  had  a  bearing  on  the 
experimental  results.  Figure  28  shows  the  effe  ct  of  three  different  waveguide  heights 
on  the  loss  characteristic  of  diode  G19-2F.  Note  that  the  "static"  or  zero  bias  insertion 
loss  is  also  affected  by  mount  geometry  and  diode  position. 
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TOP  VIEW,  SHOWING  PLACEMENT  OF  DIODE 
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SECTION  A- A  ' 7  X  taper 

(a)  Top  and  cross-section  view  (simplified) 


(lj)  Photograph 


Pig.  25.  Experimental  diode  mount. 


0080  in 


Fig.  2G.  Loss  characteristic  of  diode  G19-D2 
for  various  mounting  positions. 

In  the  preceding  discussion,  phase  shift  was  not  mentioned  because  the  nature 
of  the  loss  characteristic  provides  an  indication  of  the  phase  characteristic,  and  the 
loss  data  are  easier  to  obtain  in  the  laboratory. 

Diodes  of  the  type  yielding  the  data  of  Fig.  29a  were  the  first  to  provide  a 
significant  amount  of  controlled  phase  shift  at  140  GHz.  These  diodes  were  processed 
from  50  ohm-em  silicon  by  using  a  high  temperature  diffusion.  They  have  a  nominal 
width  of  12  mils.  While  the  phase  shift  performance  is  attractive,  the  loss  is  very 
high  even  if  the  2-1/2  dB  of  "breadboard"  diode  mount  loss  is  taken  into  account. 


The  diode  of  Fig.  29b  has  the  same  geometry  as  that  shown  in  Fig.  29a  but 
displays  a  substantially  reduced  loss  characteristic.  ALowing  for  the  breadboard 
mount  loss,  the  diede  loss  is  seen  to  be  less  than  1  dB  in  the  off”  state  (0  mA)  and 
less  than  2  dB  in  the  "on"  state  (500  mA).  This  group  of  diodes  was  processed  from 
1000  ohm-cm  silicon,  a  much  higher  resistivity  than  in  the  diode  described  above. 
The  high  resistivity  is  necessary  in  order  that  the  diode  I  layer  appear  as  a  low-loss 
dielectric  to  the  millimeter  waves  when  the  diode  is  not  biased.  Only  a  few  samples 
of  the  low-loss  type  of  diode  depicted  in  Fig.  29b  were  obtained.  Further  samples 
could  not  be  obtained  because  of  excessive  time  lag  in  material  procurement. 
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Fig.  27.  Loss  characteristic  of  diode  G19-1F  for 
various  mounting  positions. 
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Fig.  28.  Loss  characteristic  of  diode  G19-2F 
for  various  waveguide  heights. 
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PHASE  SHIFT 


LOSS 


Fig.  29.  Measured  phase  shift  and  insertion  loss  for  two 
different  PIN  diodes  at  140  GHz. 


Section  IV 


FOUR-ELEMENT  ARRAY 


In  order  to  clearly  demonstrate  the  distributed  PIN  diode  phase  shifter  in  opera¬ 
tion,  a  four-element  linear-array  antenna  and  four  4 -bit  diode  mounts  were  designed 
and  constructed.  Assembly  of  multibit  phase  shifters  and  antenna  testing  were, 
however,  precluded  by  the  substantial  insertion  loss  of  the  diodes  available  for  use 
in  the  array. 

Selection  of  a  four-element  array  was  based  on  the  need  to  limit  cost  and 
complexity.  A  straightforward  corporate  (waveguide)  feed  arrangement  of  more  than 
four  elements  in  RG138/U  waveguide  would  be  difficult  to  make  to  a  reasonable 
size,  and  would  be  too  lossy  and  too  expensive.  Figure  30  shows  the  four-element 
array  with  and  without  phase  shifters  installed.  The  waveguide  openings  on  the 
front  face  are  0.  04  by  0.  08  in  and  are  spaced  1.  84  in  apart  on  centers.  This 
geometry  will  yield  a  multilobe  radiation  pattern  with  eaeh  lobe  having  a  half-power 
be  am  width  of  about  0.  6°  . 

Figure  31  shows  the  complete  pattern  of  the  array  and  a  i  expanded  display  of  the 
region  about  the  central  lobe,  as  computed  by  standard  methods.  Although  1.84  in 
is  a  rather  substantial  element  spacing,  it  is  the  smallest  spacing  that  could  be 
obtained  without  inordinate  effort  in  the  fabrication  of  the  feed.  The  commercial 
hybrid  tees  that  were  used  \  ere  shortened  as  much  as  feasible,  and  special  short 
radius  bends  were  made. 

The  phase  shifters  are  removable  from  the  array  in  order  to  evaluate  them 
properly.  The  broadside  radiation  pattern  of  the  array  would  be  measured  without 
phase  shifters  in  the  feed  section,  as  in  Fig.  30a.  Then  the  phase  shifters  would  be 
installed  and  pattern  measurements  made  at  broadside  and  at  various  steering  angles 
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Seetion  V 


CONCLUSIONS 


The  computed  RF  characteristics  of  a  silieon  slab  of  varying  conductivity,  shown 
in  Fig.  5,  indicate  that  a  phase  shifter  with  attractive  performance  (i.  e. ,  100°  /dB  or 
better)  is  possible.  The  measured  data  on  bulk  silieon  slabs,  shown  in  Fig.  14, 
tend  to  eonfirm  the  theoretical  predictions  and  indicate  that  there  are  no  fundamental 
reasons  why  a  conductivity  modulation  phase  shifter  ean  not  be  built  that  will  offer 
attractive  performance  at  140  GHz  or  other  millimeter  wave  frequencies.  The 
measured  data  for  PIN  diodes  has  the  same  form  as  the  theoretical  data.  However, 
the  similarity  is  only  qualitative,  primarily  beeause  the  PIN  diodes  tested  did  not  fill 
the  entire  height  of  the  waveguide,  and  the  conductivity  modulation  effect  was  incom¬ 
plete  in  that  it  did  not  extend  through  the  entire  I  layer  due  to  limited  excess  carrier 
lifetime. 

The  work  described  in  this  report  has  demonstrated  the  potential  feasibility  of 
the  distributed  PIN  phase  shifter.  In  view  of  the  promise  shown  by  the  device,  it  is 
recommended  that  further  development  of  construction  techniques  (including  stacking 
of  diodes)  be  carried  out  in  order  to  build  phase  shifters  having  attractive  phase- 
shift/loss  ratios. 


Appendix  A 

COMPUTER  PROGRAMS  FOR 
PROPAGATION  ANALYSIS 
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PROGRAM  ROOT 
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PROGRAM  WAVE 


1  PROGRAM  4AVE 
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PROGRAM  WAVE  (cont. ) 
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c 

5*1 

4 

J  =  0 

54 

9LIM«0.01 

55 

M(  1)  »RJ«MR 

5*> 

H(  j  )*C J*nR 

57 

M(3)=-RJ*HR 

50 

H( 4 ) s-C J*H9 

59 

6 

J  =  J  ♦  1 

69 

R 

‘  j  I  «  o 

61 

9'J  5  1*1*4 

6? 

r  ALL  WF(9AHMA,-4(I  ),Q(I  ),WFr.) 

61 

T  F ( 0 ( 1 ) • 1  T.QLCM)  G9  Tn  3 

64 

9  0  T9  6 

65 

9 

onw.Qt  i ) 

64 

*1 1  *  I 

67 

5 

continue 

60 

JF(NI.EQ.O)  50  T9  7 

69 

iF(QiaW.LT.QLlM)  01 JM.OLOW 

79 

5AMMAbGAmMA+H{ N I 5 

7 1 

90  TH  « 

7? 

7 

TFU.E0.NIT)  GO  TO  10 

7} 

no  9  1*1*4 

74 

9 

M(i)«H(  n/io.o 

70 

90  TO  6 

76 

10 

IF  (FSTSWT.NE.O)  50  TrJ  15 

77 

FSTSWT.l 

70 

C 

79 

C 

-/RITE  HEADERS  ANO  REFFRENCF  VALUES 

89 

C 

81 

IFREO-F 

82 

WRITE(6*101)  IFREQ 

89 

TRsT/A 

84 

WR!TE(6,102>  T  R 

80 

WRITEI6* 103) 

86 

AlPHRF*  REAL  {GAMRA)fO‘,08683 

87 

BETARF»AIMAG(GAMRA)*0. 573 

88 

CALL  WF(GAMM4*CZ»W*WF0) 

89 

FREAL*  REAL  MF5 ) 

99 

FCHPX*AIRAG(  */F5) 

91 

-'RITE  (6*  104)  SIGRA»ALPHRF;«ETAPE,FREAL*ECMPx 

92 

WRITF (6* 105) 

99 

LCNT-10 

94 

C 

90 

C 

WRITE  CO^PJTFO  vHUES 

96 

c 

alpha*  re  al  < 5armm*o.08683-alpwrf 

97 

15 

98 

0ETA,AIMAG(GAMRA)*O,573-BETARF 

99 

CALL  WF(r,AHMA*CZ*WfWFG) 

ion 

FREAL*  RFAt ( WF5 ) 
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PROGRAM  WAVE  (cont. ) 


101 

rCMPX«AI'iAG(  yJF  S  > 

10? 

WR!TF(^#1q4)  SIGMA jALPHA,BFTA,FRFAL*FCMPX 

10* 

L  C  NT 

■LCNT+i 

104 

!F!LcNT,FQ.5S)  WRjTE(fe»1061 

105 

c 

1 0<s 

c 

increment  sigma 

107 

c 

10* 

I  F  ( 

ABS ( 5 I5MA-1 ,  )  .LT.O'.Ol  ) 

S  T  NIC  *  1  . 

lo^ 

I  F  ( 

ABS  (  S  I  GMA-10.  )  .LT'.O.  1  ) 

STNr.10. 

m 

T  F  ( 

ABS<SIGMA„100.  > . L T i  .  ) 

STNr-i00. 

in 

IF( 

ABS(SIGMA-1000.).LT.10' 

)  SlNf»lOOO. 

u? 

I  F  { 

ABS(SISMA-10000.  LT.IOO.  )  STNC»  10000, 

113 

lF< 

ABS(SISMA-100000.  I'.LT.i'OOO.  1  GO  TO  25 

114 

SJGMA.SIGMA+SINC 

115 

A<Z  = 

FR-CJ*SIGMA/{8.B54E-12*0MEGA ) 

u* 

CALL 

WF(GAMMA,CZ»QLOW,*IFG> 

117 

GO  Tn  4 

11* 

?1 

STOP 

119 

fno 

PROGRAM  WAVE,  SUBROUTINE  RARITZ 


1  SUBROUTINE  RAR!TZ<TI,f,AMMA".NniT> 

?  complex  gamma 

1  N  0  I  T  »  0 

4  REA0I5# IOOjEND.99)  T  I  *  BETA 

■5  1  00  P^RMAT(FA.1MX#F7.D 
*  gamma. CMPLX<  a.gjBETA) 

t  Return 

R  99  M D I T ■ 1 
9  RETURN 

19  FND 


PROGRAM  WAVE,  SUBROUTINE  WF 


1  S  JBRnUTI  ‘JE  WF(5A^MA,H,tf,WFr.) 

2  C 

3  C  COMPUTES  tp(GAMM&> |**2 

4  C 

3  C3MM0N/CMPX/AKZ#CJiRJ 

3  C3MMHN/RPL/P1# AKZS3iT,A,D 

7  C3MPLE*  ARG1#  AR52j»GAMMA»H,  A<Zj  WFGfCJiRJ 

B  REAL  PT(4)*A<ZS0>T>A;D*W 

9  A*Gl«CSQRT  (( G AMMA*H ) * ( GAMM A+H ) ♦ AK Z SO ) 

10  Ahg2«CSQRT  (  (GAMHA*H>*(GAMMA4.HI  +  AKZ*AKZsQ) 

11  WFg*ARg1*CSI'J  (ARG2*tWCCOS  (  A  R  G  2  *  T  5 

12  1  ♦ARG2*CSIM  (ARGUOI/CCOS  (  A  R  G 1  ♦  0 ) 

w*  RFALHFfi)*  REAl(WFG)*AIMAG(WFG)*AIMAG(»/FG) 
!4  RETURN 

15  FNO 
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